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Fate of baby radio galaxies: Dead or Alive ? 
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In order to reveal the long-term evolution of relativistic jets in active galactic nuclei (AGNs), we examine the dynami- 
cal evolution of variously- sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects 
(MSOs), Fanaroff-Riley type H radio galaxies (FRHs)]. By comparing the observed relation between the hot spot size and 
the linear size of radio source with a coevolution model of hot spot and cocoon, we find that the advance speed of hot 
spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRH phase). 
The deceleration is caused by the growth of the cross-sectional area of the cocoon head. Moreover, by comparing the hot 
spot speed with the sound speed of the ambient medium, we predict that only CSOs whose initial advance speed is higher 
than 0.3-0.5c can evolve into FRHs. 
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1 Introduction 

According to the standard model for radio-loud AGNs, a 
pair of relativistic jets transport away the bulk kinetic en- 
ergy from the central compact region close to a central su- 
permassive black hole (SMBH) to /h > 10 kpc scale ra- 
dio lobes, where /h is the distance from the centre of the 
galaxies (Blandford & Rees 1974). However, how young ra- 
dio loud AGNs evolve into powerful extended radio sources 
(e.g., FRn radio galaxies) is one of the primary problems in 
astrophysics. In order to clarify this issue, it is important to 
understand the nature of the small, young progenitors of FR 
Us. Young and compact radio sources such as CSOs (/h < 
Ikpc) and MSOs (/h = 1 — 10 kpc) have been discovered by 
many authors (e.g., Wikinson et al. 1994; Fanti et al. 1995; 
Readhead et al. 1996). The bright CSOs and MSOs show the 
FRn like morphology and possess hot spots (reverse shock) 
which are the signatures of supersonic expansions. The re- 
cent observations of several CSOs and MSOs have indicated 
that their age (~ 10^~^ yr) is shorter than the typical age of 
FRHs (e.g., Parma et al. 1999), which implies that CSOs 
and MSOs are possible candidates as the progenitors of FRI 
Is (e.g., Owsianik, Conway & Polatidis 1998; Taylor et al. 
2000; Giroletti et al. 2003; Polatidis & Conway 2003; Nagai 
et al. 2006). 

A number of authors have investigated the long-term 
evolution of extragalactic radio sources in different ways 
(e.g., Begelman & Cioffi 1989: hereafter BC89; Fanti et 
al. 1995; Begelman 1996; Readhead et al. 1996; Kaiser & 
Alexander 1997; Snellen, Schilizzi & van Langevelde 2000; 
Perucho & Marti 2002; Kawakatu & Kino 2006, hereafter 
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KK06). A constant advance speed of hot spots, or a con- 
stant aspect ratio of the cocoon (i.e., a self- similar evolu- 
tion) have often been assumed for a dynamical evolution of 
radio-loud AGNs for simplicity. Do these assumptions re- 
flect the actual evolution of radio sources ? In order to an- 
swer this question, it is essential to build up an appropriate 
model of radio sources without assuming the constant as- 
pect ratio and the constant advance speed of hot spots such 
as that presented by KK06. 

Recent observations suggested that the power law index 
for the evolution of hot spot size changes at the transition 
between the interstellar medium and intergalactic medium, 
i.e., /h ~ 1 — 10 kpc (Jeyakumar & Saikia 2000: here- 
after JOO; Perucho & Marti 2003: hereafter PM03). Since 
hotspots are identified as the reverse- shocked region of the 
decelerating jet, the evolution of hot spot size could reflect 
the dynamical growth of radio sources of various scales. 
However, it was hard to derive the dynamical evolution of 
radio sources from previous work (JOO and PM03) because 
of lack of spatial resolution, the observational bias and small 
sample of radio sources. In this study, we first compile sam- 
ple of CSOs, MSOs and FRIIs sources larger than in previ- 
ous works, by considering the observational bias and being 
careful about the data quality. Then, from the direct com- 
parison with KK06, we make clear a long-term evolution of 
advance speed of hot spots (Kawakatu, Nagai & Kino 2008 
for details: hereafter KNK08). 

2 Observed ths — relation 

Following the recent work of Nagai (2007), we compiled 
117 variously- sized radio galaxies and examined the physi- 
cal quantities of the hot spots. Based on these data, we focus 
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Distance from the core [kpc] 

Fig. 1 The relation of hot-spot size (rns) and hot-spot 
distance from the core (/h). Crosses with arrows indicate 
upper limit. Solid line corresponds to the best-fit for the 
sources where 10~^ kpc< /h < 1 kpc whereas broken 
line corresponds to that for the sources where 1 kpc< Ih < 
10^ kpc. Note that upper-limit data were not included in the 
fittings. 

on how the size of the hot spot changes with the distance 
from the core. Figure 1 shows the size of a hot spot, rns 
with respect to /h, together with the best linear fit on the 
log-log plane (log rns = cl log /h + We estimated rns as 
(^maj + ^min)/2 whcrc 6>maj and ^min are full width at half 
maximum (FWHM) of major and minor axes of the compo- 
nent obtained by Gaussian model fittings, respectively. The 
slope changes around ~ 1 kpc, and then the best fit value 
of a = 1.34 ± 0.24 for Ih < I kpc, while a = 0.44 ± 0.08 
for /h > 1 kpc. This tendency is consistent with PM03, but 
a large number of sources in our sample allow us to confirm 
the tendency more clearly. Concerning the possible uncer- 
tainties of the estimation of the hot-spot size and the slope 
change are discussed in KNK08 with more details. As a re- 
sult, even if we allow for possible uncertainties, the trend 
of the observed rns — relation would not be changed. It 
seems reasonable to suppose that the slope change occurs 
around 1 kpc. 



3 Coevolution model of hot spot and cocoon 

We briefly review here a dynamical evolution model of ra- 
dio sources (KK06) which traces the dynamical evolution of 
advancing hot spots and expanding cocoon (e.g., BC89; see 
also Kino & Kawakatu 2005). We consider a pair of initially 
relativistic jets propagating in an ambient medium with mat- 
ter density (pa)- Here is given by Pa(^h) oc l^^. The ba- 



sic equations of the cocoon expansion can be obtained as 
follows: 

1. The equation of motion along the jet axis, i.e., the mo- 
mentum flux of a relativistic jet is balanced by the ram 
pressure of the ambient medium spread over the effec- 
tive cross-sectional area of the cocoon head, that is, Ah, 
Lj/c = Pa^HS^h, where Lj, c and vhs are the total ki- 
netic energy of jets, the light speed and the hot spot ve- 
locity, respectively. Here we assume that Lj is constant 
in time. 

2. The equation of motion perpendicular to the jet axis, that 
is the sideways expansion velocity, Vc, which is equal 
to the shock speed driven by the overpressured cocoon 
with internal pressure (Pq), Pq = Pa^c • 

3 . The energy conservation in the cocoon, namely all of the 
kinetic energy transported by the jets is deposited as the 
cocoon's internal pressure, that is, PqVc = 2(7c — 
where Vc is the volume of the cocoon, t is the life time 
of the source and 7c = 4/3 is the specific heat ratio of 
the relativistic plasma in the cocoon. 

The cross-sectional area of the cocoon body is given by 
Ac{t) = Till (X t^, where 1^ = . Vc{t')dt' is the ra- 
dius of the cocoon body. Here tmin is the time at which the 
two-dimensional (2D) phase (Ah growth phase) starts. As- 
suming that Ah / r^s is constant in time, rns and ^hs can 
be described in terms of the length 

rns oc /h^% (1) 
^HS oc /h*^". (2) 

where S, = [X(-2 + 0.5a){a - 2) + 3a - 4]/[2X(-2 + 
0.5a) + 6] and 5'v = [2 - X{2 - 0.5a)]/[X(-2 + 0.5a) + 
3]. For a free parameter X, we can constrain the value of 
1.2 < X < 1.4, by comparison with numerical simulations 
(Scheck et al. 2002; Perucho & Marti 2007). 

4 Evolution of the advance speed of hot spots 
4.1 Constant velocity evolution ? 

In the current models of the evolution of powerful radio 
sources, the constant velocity of hot spots has often been as- 
sumed (e.g., Fanti et al. 1995; Begelman 1996; Readhead et 
al. 1996; Kaiser & Alexander 1997). But, a deceleration of 
hot spots may take place via a strong interaction with denser 
ambient gas in host galaxies since CSOs and MSOs will 
have a significant interaction with the ambient medium as 
they propagate through it (e.g., Gelderman & Whittle 1994; 
de Vries, Barthel & O'Dea 1997 ; Axon et al. 2000; Holt, 
Tadhunter, & Morganti 2008; Labiano 2008). Thus, it is still 
unclear whether vus = const is a reasonable assumption on 
the evolution of radio sources. In order to test the validity of 
the ^HS = const, model, we will derive the required mass 
density profile to explain the observed rns — diagram. 
From eq. (2), the relation of S'v ex 2 - X{2 - 0.5a) = 
is needed to realize the constant velocity of hot spots. By 
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Fig. 2 top) Growth rate of cocoon body (Ac oc t^) against 
the value of exponents of the ambient matter density profile (pa oc 
l^^). The dark grey (blue in the on-line version) region shows the 
allowed region given by the best fit value of the observed rns — W 
relation for the CSO-MSO phase (/h < 1 kpc), i.e., = 1.34 zb 
0.24. The thick dotted line (red thick dotted line in the on-line 
version) corresponds to the allowed region obtained for the MSO- 
FRH phase (/h > 1 kpc), including the fitting errors, i.e., — 
0.44 zb 0.08. The black solid line shows ^;HS=const (5'v = 0). The 
two black dashed lines are relations for the evolution of hot spot 
velocity ^;hs oc . The filled circles show the solutions for our 
model. We can constrain the power-law index a(/h < 1 kpc) as 
< a < 0.5. The filled triangles are solutions for the ^;hs = 
const., bottom) The predicted ambient mass density profiles are 
shown in the bottom panel. 



eliminating X in eq. (1), one finds rns l^^^ . Compar- 
ing this equation with the observed rns — relation, the 
constant velocity model requires that the slope of the inner 
part of the density profile (< 1 kpc) must be steeper than 
that of the outer part (> 1 kpc), i.e., a ~ 2.6 (/h < 1 kpc) 
and a ~ 0.8 {Ih > 1 kpc), as seen in Fig. 2 bottom (filled 
triangles and also the schematic picture). However, such a 
density profile of ambient matter is unrealistic because of 
the slope of the density profile in many clusters of galax- 
ies and groups of galaxies (/h > Ikpc) where a ^ 1.5 
(e.g., Trinchieri, Fabbiano & Canizares 1986). Thus, the 
^HS = const model can be rejected. 

4.2 Deceleration and acceleration of hot spots velocity 

We show an evolutionary track of radio sources that is con- 
sistent with the observed evolution of rns (see Fig. 1). In 



groups of galaxies and clusters of galaxies (/h > 1 kpc), 
it is well established that the slope of the ambient matter's 
density profile is a ^ 1.5 (e.g., Trinchieri et al.l986). We 
assume here a = 1.5 for /h > 1 kpc. Comparing the KK06 
model with the observed rns — relation, we can deter- 
mine (i) the evolution of the advance speed of hot spots and 
(ii) the slope of mass density distribution for /h < 1 kpc. 
The predicted evolution of hot spots shows that the advance 
speed of the spots and lobes shows the deceleration phase 
(CSO-MSO phase) and the acceleration phase (MSO-FRD 
phase) as follows (see filled circles showing allowed solu- 
tions in Fig. 2): In the CSO-MSO phase (/h < 1 kpc), the 
hot spot decelerates as 



^HS OC l^^. 



(3) 



In the MSO-FRU phase (/h > 1 kpc), the hot spot slightly 
accelerates as 



^HS OC 



0.3 



(4) 



For the CSO-MSO phase (/h < 1 kpc), a flatter density 
profile (0 < a < 0.5) is predicted in order to satisfy both the 
observed rns — relation and the constraints from numer- 
ical simulations. The mass density profile is quite similar 
to a King-profile, as indicated by X-ray observations of el- 
liptical galaxies (e.g., Trinchieri et al. 1986). The predicted 
mass density distribution implies that a significant interac- 
tion of the AGN jets with the ambient medium can take 
place during the CSO-MSO phase, as compared with the 
MSO-FRn phase. Then, the interaction between the jets and 
the ambient medium is stronger in the CSO-MSO phase. 
This leads to the larger velocity of sideways expansion in 
order to maintain the energy conservation in the cocoon. 
Thus, Ah (or rns) grows faster in the early phase of evo- 
lution. From the equation of motion along the jet axis, it 
is found that the advance speed of hot spots (vus) is de- 
termined by the linear density of the effective working sur- 
face, Pa^h, i.e., 'i^HS ^ (Pa^h)""^. When pa^h incrcascs 
with /h in the CSO-MSO phase (pa^h oc l^), the hotspot 
velocity decelerates. On the other hand, the advance speed 
increases if Pa^h decreases with /h in the MSO-FRII phase 
(Pa^h oc ^■^). Summing up, the deceleration and accel- 
eration of hot spot velocity is caused by the deviation from 
the balance between the deceleration effect via the growth 
of cocoon head and the acceleration effect due to the de- 
crease of the ambient mass density. 

5 Prediction of fate of CSOs: Dead or Alive ? 

We further investigate which kind of CSOs can evolve into 
FRn sources. For this aim, we compare the predicted evolu- 
tion of ^HS with the sound velocity of the ambient medium, 
Cs, since the supesonic jets can maintain at /h > 10 kpc 
only when vns > Cg. For the slope of ambient matter den- 
sity, we assume a (/h < 1 kpc) = and a{l\i > 1 kpc)=1.5 
(see §4.2). Correspondingly, the behavior of ^hs can be de- 
termined as ^HS oc l^^ for /h < Ikpc and vus oc l^-^ 
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Fig. 3 Our predictions of evolution of vhs, i.e., ^hs oc 
l^^ for /h < Ikpc and vyls oc /{|-^ for /h > 1 kpc. 
The black solid lines denote the evolution of vus for 
'^Hs(^h,2D) = O.Olc, O.lc, and 0.5 c at /h, 2D = 5pc. The 
light grey shaded region represents the evolutionary path 
from CSOs into FRIIs. The black shaded region shows 
the range of sound velocity of the ambient medium, i.e., 
7xl0-^c < Cs < 1.4xlO-^c(5xlO^K < Tg < 2xlO^K) 
where Tg is the temperature of the ambient medium. 



for /h > 1 kpc (see eqs. (3) and (4)). The hot ambient- 
gas temperature, Tg is measured to be in the range of Tg = 

5 X lO^K - 2 X 10^ K, i.e., c, = {bkTj'^m^Y^^ ^7 x 
IQ-^c - 1.4 X IQ-^c (e.g., Trinchieri et al. 1986), k is the 
Boltzman constant and rrip is the proton mass. 

Figure 3 shows our predictions of the evolution of hot 
spot velocity for three initial advance speeds in the ambient 
medium with VHs(^h,2D) = O.Olc, O.lcandO. 5c. Here we 
suppose /h,2D = 5pc where /h,2D = ^Hs(^h,2D)^min is the 
distance from the core at which the 2-D phase starts. As seen 
in Fig. 3, we find that CSOs can evolve into FRII sources, 
passing through MSOs for any /h when VHs(^h,2D) is larger 
than 0.3-0. 5c because of VHs(^h) > Cg. On the other hand, 
when VHs(^h,2D) is less than 0.3-0. 5c, vhs is comparable 
to the sound velocity during the CSO-MSO phase (/h < 
1 kpc). Thus, the CSOs can evolve into FRII sources, pass- 
ing through distorted MSOs. Such distorted MSOs might 
correspond to low power compact radio sources of sizes 
of ~ kpc (e.g., Kunert-Bajraszewska et al. 2005; Giroletti 
2007). 

6 Summary 

We have investigated the relation between CSOs, MSOs and 
FRIIs, by comparing the coevolution model of hot spots and 
a cocoon (KK06) with the observed rns — relation reflect- 
ing the dynamical evolution of radio sources. We find that 
the advance speed of hot spots and lobes strongly deceler- 
ate when the jets pass through the ambient medium in host 
galaxies (i.e., the CSO-MSO phase), while the jets accel- 
erate outside host galaxies (i.e, the MSO-FRE phase). The 



reason of deceleration is the growth of the cross-sectional 
area of cocoon head. The predicted deceleration of hot spots 
seems to be consistent with the recent observation show that 
the outflow velocity of the ionized gas around radio lobes 
may decelerate as /h increaes (Labiano 2008). Furthermore, 
by comparing the hot spot speed with the sound speed of the 
ambient medium, we predict that only CSOs whose initial 
advance speed is higher than 0.3 -0.5c can evolve into FRI 
Is. This indicates that the origin of the FRI/FRE dichotomy 
could be related to the difference in the initial advance speed 
(Kawakatu, Kino & Nagai 2008 in preparation). 
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